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Abstract. In this paper we are concerned with various graph invari- 
ants (girth, diameter, expansion constants, eigenvalues of the Laplacian, 
tree number) and their analogs for weighted graphs - weighing the graph 
changes a combinatorial problem to one in analysis. We study both 
weighted and unweighted graphs which are extremal for these invari- 
ants. In the unweighted case we concentrate on finding extrema among 
all (usually) regular graphs with the same number of vertices; we also 
study the relationships between such graphs. 



1. Introduction 

In this paper we are concerned with various invariants of graphs related 
to distance (girth, diameter), expansion properties of the graph (Cheeger 
constant) and spectrum of the Laplacian (the first eigenvalue, the tree num- 
ber). 

The behavior of many of these invariants for random graphs has been 
extensively studied ( []Bo| ). For regular graphs, the results about weak 
convergence of graph spectra to the spectrum of the infinite regular tree 
( flKest| , |Mck81| ) allow to get many results on the asymptotic behavior of 



spectral invariants (as the number of vertices increases). 

Often when certain bounds are established for the rate of growth of some 
invariant, one is interested in studying families graphs with (asymptotically) 
the best possible behavior with respect to a given invariant (two well-known 
examples are Ramanujan graphs and cages). 

A related problem is to study the graphs which maximize (or minimize) 
a given invariant among all simple connected graphs with n vertices and 
m edges (we denote the set of such graphs by G n ,m), or among fe-regular 
connected simple graphs with n vertices (we denote the set of such graphs 
by <-fn,fc)- The sets Q nim and X n ^ are finite, so maximizing a given invariant 
is a discrete optimization problem. 

On the other hand, Q n ^ m (respectively, X n ,k) is naturally embedded in 
the complex Gnm °f an weighted graphs with n vertices and m edges (re- 
spectively, the complex X™ k of weighted £>regular graphs with n vertices); 
the edge weights are normalized so that their sum is constant. One can 
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study the extremal problems for weighted analogs of graph invariants (see 
e.g. [jFiefl ). A natural question then is which unweighted simple graphs in 
Qn,m or X n k (all weights are equal to 1) are local extrema in Q™ m or X™ k 
for a given invariant. In section ^ we study this question for various graph 
invariants. 

In section [3| we discuss some heuristics about the structure of graphs which 



have many spanning trees. In section 4.1 we give an upper bound on the 



rate of growth of the number of spanning trees for 6-regular simple graphs 



on a fixed surface. In section 4.2 we prove a result about the asymptotic 
properties of certain lattices, associated to graphs by Bacher, de la Harpe 
and Nagnibeda (cf. pHNj ), for sequences of regular graphs. Finally in sec- 
tion H we explain how the first 8 cubic cages appear naturally in a sequence 
of graphs produced from by "greedy" edge insertion with respect to the 
number of spanning trees. 

In the remainder of this section we review the definitions and some results 
about the asymptotic behavior of various graph invariants. 

1.1. Definitions and notation. The adjacency matrix A of a graph G € 
G m ,n is a square matrix of size n where Aij is the number of edges joining 
the vertices u, and vj. We always consider loopless graphs, so An = 0; for 
simple graphs Aij is equal to or 1. 

The nearest neighbor Laplacian L acts on functions on the set V(G) of 
the vertices of G: given / : V ->■ R, L(f)(v) = Y, w ~v(f( v ) ~ f( w ))- Let 
L = {Lij} be the matrix of L(G) of a (not necessarily simple) graph G; 
then La is the degree of the vertex v i, and —Lij is the number of the edges 
joining V\ and Vj ^ Vi (equal to or 1 for simple graphs). For fc-regular 
graphs L = k ■ Id — A. 

Let A be the adjacency matrix of a graph G G X n ,k, and let its spectrum 
(in the decreasing order) be given by 

k = Ai > A 2 > A 3 > . . . > A n > (-k) (1) 

(the constant function is the eigenfunction with the eigenvalue k). The last 
eigenvalue A n = —k if and only if G is bipartite. The spectrum of L(G) is 
Hj = k — Xj+i,j = 0, 1, . . . , n — 1; it satisfies = no < Hi < . . . < fJ> n ~i < 
minjra, 2k}. 



1.2. Spectral properties of random graphs. Alon and Boppana ( [ A186 , 
LPS|| ) proved that 



lim inf X 2 (G) > 2 ^/k^l (2) 



(cf. |Nil] , Bole] for more precise results); /c-regular graphs all of whose Xj ^ 
zkk lie inside [— 2yk — 1, 2\/k — 1] are called Ramanujan graphs (because 
of the relation to Ramanujan conjectures, cf. [ |Sar9C , Lub| ]). The first such 



graphs were constructed in [ LPS , Mar8£ ] ; many other constructions followed 
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1.3. Complexity of a graph. An important invariant of G is the number 
t{G) of spanning trees of G; it is sometimes called the complexity of G. By 
Kirckhoff's theorem ( flKirfl ), 

n t(G) = • . . . • 

and t(G) is equal to the determinant of any cofactor of the matrix L and 
(it is some times called the determinant of Laplacian). McKay proved in 
[|Mck83| that 

lim ( sup (r(G)) 1 /^ = a k = f - (3) 

and that a typical graph G in X n ± satisfies 

(r(G)) 1/n ~ a k 

The above formula also holds for some Ramanujan graphs ( |Sar|| ). 



Alon in [ A19Q | proved some results concerning 

* " Jfe Ut (T(G))1/ ") (4) 

Namely, he showed that for k > 3, S k > y/2 and that 5(k)/k = 1 + o(l) as 



k — ► oo. Valdes in | VaJ j determines the 2-connected cubic graphs with the 
minimum number of spanning trees. It follows from her results that if the 
limit in (0) is taken over all 2-connected cubic graphs, the answer is 

2 3/4 



(this was also proved by A. Kostochka in | Kost |). 



1.4. Extremals for girth. We begin with girth. It is easy to show that a 
fc-regular graph (k > 3) with girth g > 3 has at least fo(k,g) vertices where 

j l + k^k-l)^' 2 -\}/{k-2), godd 
Uk,9) = { 2p-l)^-l]/(fc-2), 9 even (5) 

The fe-regular graphs with girth g and fo(k,g) vertices are called Moore 
graphs (sometimes they are also called generalized polygons for even g > 4) ; 
these graphs are unique if they exist. For g = 3 those are complete graphs 
K k+ i, for g = 4 those are complete bipartite graphs K k ^. It is proved in 
[ |Big93| , Ch. 23] that for g > 4 Moore graphs only exist if g = 5 and k = 3 
(Petersen graph), k = 7 (Hoffmann-Singleton graph), and (possibly) k = 57; 
or if g = 6,8 or 12. /c-regular Moore graphs of girth 6 exist if and only if a 



projective plane with k points on the line does [ Sin 



The /c-regular graphs of girth g which have the smallest possible number 



of vertices are called (k,g) cages (such graphs always exist, cf. [HS, Ch. 
6]), and the number of their vertices is denoted by f(k,g). For k > 3 and 
g such that the corresponding Moore graphs exist, those graphs are (k,g) 
cages f(k,g) = fo(k,g); if Moore graphs don't exist, f( k,g) > fo(k,g). 
Very few cages other than Moore graphs are known (see |Won| 1, [HS, Ch. 



6] and [Big9S, Ch. 23] for surveys of known results). For cubic graphs, 
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f(3,g) is known for 3 < g < 12 ( flBMSj ). For k > 4, the f(k, 7) is already 
unknown. The best upper bounds for f(k,g) come from various infinite 
series of symmetric graphs (cf. |LPS| , [Mar88| , pVIorj , Chi , [BH| , |LUW| ). It is 
known that 

1/2 < limmf l0Sfc - l/(fc ' g) < limsu P l0g *- l/(fc ' g) < 3/4 (6) 
g g^oo g 

For all known infinite families of regular graphs, 

limsup^og^^! f(k,g))/g = 3/4 

g^oo 

it is sometimes conjectured that the upper bound in (^) is actually an as- 
ymptotic result. 



1.5. Extremals for complexity. We next turn to the number of span- 
ning trees. It is easy to see that the complete graph on n vertices has the 
maximum number of trees among all multigraphs with the same number of 
vertices and edges. In a series of papers (see e.g. |CK|, |Kel74|, |Kel76|, |Kel80|, 



Kcl96, Kcl97] 1) Kelmans characterized the graphs, obtained by removing a 
small number of edges form the complete graph, which have the maximum 
number of spanning trees among all the graphs with the same number of 
vertices and edges. Cheng ([ phe |) proved that the complete multipartite 
regular graphs also have this property; see also |Con| 



We notice that -Kfc+i and K^^ are, respectively, the unique (k, 3) and 
(k, 4) cages. One can also check (see [ |Va| ) that the Petersen graph (which 
is the (3, 5) cage) has the most trees (2000) among all cubic simple graphs 
with 10 vertices. We conjecture 

Conjecture 1.1. Moore graphs have the maximum number of spanning 
trees among all simple graphs with the same number of vertices and edges. 

Some evidence supporting this conjecture is discussed in section ||. 

Acknowledgements. The authors would like to thank N. Alon, P. Sar- 
nak, P. Seymour and R. Wilson for helpful conversations, and A. Kelmans 
for explaining many results in extremal graph theory during his visit to 
Caltech in the spring of 1997, as well as for suggesting many improvements 
to the manuscript. The first author would like to thank the Institute for 
Advanced Study at Princeton for their hospitality in the fall of 1997 and in 
the spring of 1998. 



2. Local extrema for weighing the edges 

In this chapter we shall discuss several extremal problems for weighted 
graphs. We will be concerned with finding extrema for various functionals 
on the weighted edge polytope P (which is actually a simplex) defined below. 
Let G be a simple connected graph, and let E(G) = {e%, ■ ■ ■ , e m } be the set 
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of its edges. Then 

f m \ 

P(G) : = < (xi, . . . ,x m ) : Xj >0, J^Xj = m > (7) 

where the (nonnegative) real number Xj is called the weight of e^; the cor- 
responding graph is denoted by G(xi, . . . , x m ). The graph G(l, 1, . . . , 1) is 
called an unweighted graph. We introduce the following notation: 

w = (1,... ,1). 
The adjacency matrix of G(x\, . . . , x m ) is defined by 



i = j or Vi 9^ Vj 
eu = {viVj) 



The degree deg(i>j) of a vertex Vi is defined by 

deg(^) = °°k 

j:(ViVj)=e k eE(G) 



The Laplacian is defined by 



Lij 



deg(vi), i = j 

0, Vi ^ Vj (8) 

-Xk, e k = (viVj) 



The notation for the eigenvalues of A and L will be the same as for the 
unweighted graphs. The proof of Kirchhoff 's theorem extends to weighted 
graphs ( [ |Chu | ) ; we shall continue to call the determinant r(G) of any cofactor 



of the (weighted) tree number of G. 

One can regard deg(uj) as a combinatorial analog of a curvature at «j. Mo- 
tivated by an analogy with metrics of constant curvature, we define another 
polytope P{G) C P(G) by 

P(G) := {(xi,... ,x m ): xj>0, deg^) = d> Vi} (9) 

where di = #{j / i : Vi ~ Vj} is the degree of Uj in the unweighted graph 
G. 

We will consider the problem of maximizing or minimizing various func- 
tionals over P and P. Some of the functionals (girth and diameter) will be 
linear in Xj-s, and so we will be solving a problem in linear programming 
( |Sch| , Ch. 7]). Other functionals (expansion constants) will not be linear in 
Xj-s on P, but will become linear on P. 

We shall be often concerned with the following question: 

Question: When is u = (1, 1, . . . , 1) an extremum (in P or P)? 

In other words, when is the unweighted graph an extremum for a given 
functional? In some cases we are able answer the above question completely. 
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2.1. Girth and diameter. We turn to ths study of the girth and the 
diameter. We make several definitions: The length of a path in a weighted 
graph G is the sum of the weights of the edges in the path. The girth 7(G) 
is the length of the shortest cycle in G. The distance between two vertices of 
G is the length of the shortest path connecting them; the diameter diam(G) 
is equal to the largest such distance. 

We are interested in the points x £ P(G) where 7(3;) takes its maximal 
value and where diam(x) takes its minimal value on P. Denote by C\, . . . ,C S 
all the cycles of G. Then finding the maximum of the girth leads to solving 
the following problem in linear programming: 



We assume that any edge on a cycle is visited at most once, so there are 
finitely many cycles; we make a similar assumption about the paths in the 
next definition. 

Similarly, for two vertices u, v of G denote by {iti(u, v ), . . . , n a (u, v)} the 
set of all paths in G connecting u and v (where a = a(u,v)). Finding the 
minimum of the diameter leads to the following problem: 



We now fix some x £ P(G) and consider the corresponding graph G = 
G(x). Two vertices u, v of G are called diametrically opposite if dist(«, v) = 
diam(G). Any path connecting two diametrically opposite points which has 
length diam(G) is called a meridian. Also, any shortest cycle is called a 
systole. 

To study how girth changes, we remark that for small weight changes 
the new systoles will come from among the weighted old systoles. Denote 
the systoles of G by Si, S2, ■ ■ ■ ,S r , and form the (m-by-r) edge-cycle inci- 
dence matrix Sij. Turning to the diameter, we remark that for small weight 
changes the new meridians will come from among the weighted old meridi- 
ans. Denote the meridians of G by Mi, Mi, . . . , Mi and form the (m-by-i) 
edge-meridian incidence matrix Mjj. Denote by S (respectively, M) the 
cone formed by all linear combinations of column vectors of Sij (respec- 
tively, Mij) with nonnegative coefficients. Clearly, both C and M lie inside 
the first quadrant. 

Theorem 2.1. The graph G is a local maximum for girth (respectively, a 
local minimum for diameter) if and only if to lies in S (respectively, M. ). 

Proof. Let G(x°) be the graph in question, and let x 1 be another point 
in P(G). Denote by y the vector x 1 — x° connecting the points in P. The 
graph G(x°) is a local maximum for the girth if and only if the inequality 





mm max mm 

xeP I u,veG l<j<a(u,v) 




Sy > 
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has no solutions for small enough y. By Farkas' lemma ( |Sch| , Ch. 7]), this 
condition is equivalent to the condition u> G S. 

Similarly, the graph G(x°) is a local maximum for the girth if and only if 
the inequality 

My < 

has no solutions for small enough y. Again by Farkas' lemma, this is equiv- 
alent to the condition u £ Ai. 

□ 

An obvious sufficient condition for the assumptions of the Theorem |2.1| 
to hold is for every edge to be contained in the same number of systoles 
(diameters), so that edge-transitive graphs are always satisfy this assump- 
tion. Another large class of graphs where every edge is contained in the same 
number of systoles consists of graphs arising from 4- vertex-connected proper 
triangulations of surfaces (every edge is contained in exactly two triangles). 

On the other hand, the assumptions of Theorem are clearly not sat- 
isfied when there is an edge in G which is not contained in any systole 
(diameter) . 

2.2. Expansion constants. We now turn to Cheeger and expansion con- 
stants. Let G be an unweighted graph. Given a partition of its vertex set 
G = A U B into disjoint sets A and B, define E(A, B) to be the set of all 
edges with one endpoint in A and another in B. Let \E(A, B)\ be the sum of 
the weights of the edges in E(A, -B)|] and let \A\ denote the sum of degrees 
of all vertices in A (also called the volume of A)^ The expansion constant 
c(G) is defined by 

c(0) _ mi „ H (10) 

V ; G=AuB \A\ -\B\ v ' 

The Cheeger constant h(G) is defined by 

h(G) = min 1^,?' (11) 
v ; g=aub min(|^4|,|B|) 



To prove an analog of Theorem 2.1 for the above constants, we find it 



convenient to maximize Cheeger and expansion constants over P{G) instead 
of P{G). The reason is that for A C V{G) the volume \A\ does not depend 
on x E P(G) (since the degrees are preserved), and so c(G) and h(G) become 
linear functionals on P(G). 

Let Nij denote the m-by-n edge-vertex incidence matrix. The dimension 
of P is equal torn - rank(JVy). We denote by N the cone spanned by linear 
combination of the column vectors of (Nij) with nonnegative coefficients.^] 



1 It is equal to the number of the edges in E(A,B) for unweighted graphs. 

2 Often |j4| denotes the number of vertices in A, which for regular graphs is proportional 
to our definition; our definition allows a straightforward extension to weighted graphs. 

3 Here TV stands for normalization conditions; in general, the more conditions we impose, 
the larger the cone TV is; vector uj_ should always be in A/". 
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To state the next theorem, we need to define several matrices. Given 
x G P, let G(x) be the corresponding graph. An edge cut H of G(x) is called 
an expansion edge cut (respectively, a Cheeger edge cut) if it is one of the 
cuts for which the expansion constant (respectively, the Cheeger constant) 
of G is attained. Let H\, H2, ■ ■ ■ , H r (respectively, H , H 2 , ... , H l ) be the 
expansion (respectively, Cheeger) edge cuts. 

Denote by A{H) and B(H) the disjoint subsets into which H separates 
V(G). We now define two matrices: an m-by-r expansion matrix (Ey) and 
an m-by-i Cheeger matrix (Cij) by 

= f 1/(1^)1-15(^)1), eteHj 
%3 1 0, otherwise. 

and by 

c .. = I l/mind^)!, 1^(^)1), e^GiF 
JJ I 0, otherwise. 

Let £ (respectively, C) denote the cone formed by linear combinations of 
column vectors of Eij (respectively, Cij). 
We are now ready to formulate and prove 

Theorem 2.2. A graph G is a local maximum for the expansion constant 
(respectively, the Cheeger constant) for degree-preserving edge valuations if 
and only if £ n N ( respectively, C Pi N) contains a nonzero vector. 



Proof : The proof is similar to that of Theorem 2A . □ 
We remark that edge-transitive graphs satisfy both assumptions of The- 
orem |2l| (for such graphs u G (£ fl Af) and u G (C n AA)). Also, as in the 
case of girth and diameter the assumptions are not satisfied if some edge is 
not contained in any Cheeger (respectively, expansion) edge cut. 

2.3. The tree number. By the weighted version of Kirckhoff's theorem 

( M ) 

r(G) = £ II *i ( 14 ) 

TeT(G) ej&T 

where the sum is taken over the set T(G) of the spanning trees of G. 

Finding the maximum of t(G) on P becomes a Lagrange multiplier prob- 
lem with the constraints given by Eq. (0). The condition for x G P(G) to 
be a critical point for t(G) is 

dr{G) dr{G) dr(G) 



dxi 8x2 dx m 

The partial derivatives above are given by 

dr(G) x - -pr 

3 TeT(G) k^j;e k eT 

The ratio tJt{G) is called the effective resistance of e,-. 



(15) 
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We have thus proved: 

Proposition 2.3. The graph G(x) is a critical point for t(G) if and only 
if the effective resistances of all edges are the same. 



If an unweighted graph satisfies the assumptions of Proposition then 
every edge of this graph is contained in the same number of spanning trees. 
Such graphs were studied by Godsil in | Goc| ; he calls these graphs equiar- 



boreal. Obviously, all edge-transitive graphs (the automorphism group acts 
transitively on the edges) are equiarboreal.^ Godsil gives several more suffi- 
cient conditions for a graph to be equiarboreal; in particular, any distance- 
regular graph and any color class in an association scheme is equiarbo- 
real (the least restrictive condition Godsil gives is for a graph to be 1- 
homogeneous). By an easy counting argument one can show that for an 
unweighted equiarboreal graph 

Tx = T 2 = ... = T (G)-(n-l)/m (16) 



(this is actually the result of Foster, cf. [ |Fos| ) so the necessary condition for 
a graph to be equiarboreal is 

m | (n- 1)t(G). 

We remark that the graphs which have the most spanning trees among the 
regular graphs with the same number of vertices are not necessarily equiar- 
boreal, and vice versa. For example, the 8-vertex Mobius wheel (cf. section 
||) which has the most spanning trees among the 8-vertex cubic graphs is 
not equiarboreal (cf. also |Val| |), while the cube (which is certainly edge- 
transitive, hence equiarboreal) has the second biggest number of spanning 
trees among the 8-vertex cubic graphs. 

We next study the function t(G) on P(G). Surprisingly, it turns out to 
be concave, so there is at most one critical point, and if there is one it has to 
be a global maximum. The eigenspace Eq corresponding to the eigenvalue 
/io = of L on G(x) is the same for all x in the interior of P (and consists 
of constant functions). Let LTo be the operator of projection onto Eq; its 
matrix is given by (1/ra) J, where J is the matrix of 1-s. The eigenvalues of 

M = L + n 

are {1, /xi, . . . , fi n }, hence 

u 2 t(G) = det(J + L). (17) 



This result is due to Temperley ([ Tem| ]). Thus it suffices to study the con- 



cavity of the positive-definite operator M . 
Theorem 2.4. The function t(G) is concave on (the interior of) P. 
The theorem immediately implies 



4 See [Bon | 



for examples of edge-transitive graphs which are not vertex-transitive. 
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Corollary 2.5. If z is a critical point of t{G) in the interior of P, then it 
is the global maximum for t(G). In particular, the unweighted graph is the 
global maximum for equiarboreal graphs. 



Proof of Theorem 2.4. It suffices to prove that the determinant is a 
concave function on the set of positive definite symmetric matrices. Let Q 
be such a matrix, and let 

Q(t) = Q + tB, t g R 
be a line of symmetric matrices through Q. Then 

= tr( B Q->), 

and 

d2lQgd ;; (9(t)) = -tr W -^Q-). 
at 1 

It suffices to show that the last trace is strictly positive. The matrix R = 
Q^ 1 is positive definite. Conjugate R to a diagonal matrix with eigenvalues 
w = (u>i, . . . , w n ). Then 

tr(BRBR) = wSw* 

where Sij = RijRij is the Hadamard product of R with itself. By Schur 
product theorem ( |HJ| , Theorem 7.5.3]), S is positive definite, which finishes 
the proof. □ 



We remark that Theorem 2.4 can be considered as an analogue of the 
convexity of the determinant of Laplacian in a conformal class of metrics on 
surfaces of genus greater than or eqaul to one, proved by Osgood, Phillips 



and Sarnak in [OPS|. It would be interesting to characterize graphs G where 



t(G) achieves its maximum in the interior of P(G). 

2.4. Eigenvalues of the Laplacian. We introduce some notation which 
will be used in studying the Taylor expansion of the eigenvalues. 
Given two points x , x 1 £ P(G), let y = x 1 — x° and let 

x(t) = x° + ty, < t < 1 (18) 

be the segment (call it 7r(i)) connecting the two points in P. The normal- 
ization condition reads 

ycJ = (19) 
The basis of all y-s satisfying (|l9|) is given by vectors 

hj = §i ~ % (20) 

for i ^ j where gj-s are the standard basis vectors. 

We next turn to studying how the eigenvalues of the Laplacian depend on 
t. We look at one parameter family L(i) of Laplacians on G corresponding 
to the points on 7r(t). The Laplacian is given by the matrix 

L(t) = L + tB 
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By standard results in perturbation theory ( |Kat| ) , there exists an orthonor- 
mal basis of the eigenvectors of Q such that the the eigenvalues and the 
eigenvectors admit Taylor expansions in t. 

Let fj, = fi(0) > be an eigenvalue of the multiplicity p > 1. Let be 
the corresponding (p-dimensional) eigenspace, and let P^ denote the corre- 
sponding projection. The condition for the graph G to be a critical point 
for fi is that the matrix P^BP^ be indefinite for all B. We now compute 
the eigenvalues of P^BP^ for the edge weightings given by (^) with b as in 
(f20|). We distinguish between two cases: 

a) edges e± and e2 are disjoint. 

b) edges e\ and e2 share an endpoint 

We first consider the case a). Let fx, f2, ■ ■ ■ , fp ■ V — > R be a basis of 
eigenvectors of E^. For convenience, we denote fj(v\) — fj{v2) by Zj and 
fj( v 3) — fj( v 4) by Wj. The eigenvalues of P^BP^ are given by the eigenvalues 
of the linear form F(a%, . . . , a p ) defined by 

F( ai ,...,a p ) = {B{ ai h + ... + a p f p )J 3 ) (21) 

It is easy to see that the matrix of F = Fy is given by 

Fij — Z^Zj WiWj 

The eigenvalue fi is critical if F is indefinite for any choice of B. 

We remark that F is a difference of two rank one matrices, F = z <g> z — 
w<g>w where z = (z% , . . . ,z p ) and w = {w\ w p ) and so has rank at most 
two. Computing the coefficient C p _2 of t p ~ 2 in the characteristic polynomial 
of F, we find that 

Cp-2 = ~ ( ZiW i ~ z i w i) 2 - 

l<i<j'<p 

Therefore, C7 p _2 < and is strictly negative unless x and y are proportional 
or one of them is zero. Now, if C p -2 < then F is indefinite (since it has 
at most two nonzero eigenvalues). If C p _2 = (and F has rank one) and if 
H is critical then F = 0. 

Consider first the case of a simple eigenvalue (p = 1). Then [i is critical if 
F = Fu = for every choice of B and therefore z\ — w\ = or z\ = ±w\. 
Recalling the definition of z\ and w±, we see that for the eigenvector fi = f 

\f(vi)-f(v2)\ = \f(v 3 )-f(v 4 )\ (22) 

The following proposition is proved by repeating the argument for all pairs 
of nonadjacent edges; if some edge e of G is adjacent to every other edge, we 
prove an analogue of (|22|) for the case b) of adjacent edges (it was proved 
before for the case a) of non-adjacent edges). 

Proposition 2.6. A simple eigenvalue fj, of the Laplacian on G is critical 
if and only if the corresponding eigenvector f : V — > R satisfies 



| /(it) - f(v)\ = const 



(23) 
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We now study the graphs which admit an eigenvector / : V — * R satisfy- 
ing ( |23| ) for some c > 0. If c = then / is a multiple of a constant vector and 
so has eigenvalue zero which is a contradiction. If c 7^ then it is easy to see 
that the graph G cannot have odd cycles and hence is bipartite. Namely, let 
uiu 2 ...«; bea cycle. Then (putting u t = u ) £i=i - /(ui-i)) = 0. 

But each term in the sum is equal to ±c, and since the number of terms in 
the sum is odd, they cannot add up to 0. 

We now want to study the unweighted /c-regular graphs which have an 
eigenvector (corresponding to an eigenvalue u > 0) satisfying ( |23| ) (without 
necessarily assuming that fi is simple). We shall rescale the eigenvector so 
that c = 1 in (^). From (^) and ( |23| ) it follows that for each vertex u the 
expression u • f(u) can only take one of the values k, k — 2, k — 4, . . . , —k + 
2,—k. Consider first the vertex uq where f(u) takes its maximal value a (by 
changing the sign if necessary we can assume that a > 0). It follows that / 
takes value a — 1 on all the neighbors of uq, hence 

an = k 

Next, consider any neighbor u\ of u. The value of / at any neighbor of 
u\ can be either a (let there be r\ > 1 such neighbors; no is one of them); 
or a — 2 (it follows that there are k — r\ such neighbors). From (||) it follows 
that 

fi(a — 1) = k — 2ri 
It follows from the last two formulas that 

H = 2n (24) 

where n > 1 is a positive integer. If r\ = k, then u = 2k is the largest 
eigenvalue of L. 

We next define the level of a vertex u to be equal to j if /(«) = a — j; we 
denote the set of all vertices of G at level j by Gj. It is easy to see that if 
u G has rj neighbors where / takes value a — j + 1 then 

It follows that rj is the same for all u £ Gj. Using (|24] ) we see that 

n • j = re- 
consider now a "local minimum" u £ Gn- Then = fc, and we see that 

n I (25) 

Let denote the number of vertices in Gj. Counting the vertices con- 
necting Gj-and in two different ways, we see that for all < j < N — 1, 

nj(k-rj) = n j+1 r j+ i 

Consider the case r\ = X,jx = 2. It follows from the previous calculations 
that rj = j and that N = k. Accordingly, rij = no( •) and 

|G| = 2 fc n (26) 
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We next describe a class of graphs admitting an eigenvector of L with u = 2 
satisfying (p3|). 

An obvious example of such a /c-regular graph is the fc-cube, and any such 
graph has the same number of vertices as a disjoint union of no cubes by 
(26). Start now with such a union, choose the partition of the vertices of 



each cube into "levels" and take two edges u\u 2 and U3U4 in two different 
cubes such that ui, U3 are both in level j while u%, U4 are both in level j + 1. 
If we perform an edge switch 

(uiu 2 ), (u 3 U4) -> (u 1 u i ),(u 3 u 2 ) 

then the number of the connected components of our graph will decrease 
while the eigenvector / will remain an eigenvector with the same eigenvalue. 

Performing sequences of edge switches as described above, we obtain ex- 
amples of connected graphs satisfying (|2^) and (p6| ) for any uq. Conversely, 
it is easy to show that starting from a graph satisfying ( f23| ) and (|26| ) and 
having chosen a partition of its vertices into levels one can obtain uq disjoint 
fc-cubes by performing a sequence of edge switches as above. 

We now want to consider the case when fx = /ii is the lowest eigenvalue 
of the Laplacian. The first remark is that then necessarily /x < k, and u = k 
only if G = Next, we want to consider "small" k for which k — 2\/k — 1 

(the "Ramanujan bound") is less than 2 (this happens for 3 < k < 6). It 
then follows from the results of Alon ([ [Nil| ) that the diameter of G (and 
hence the number of vertices in G) is bounded above. 

Proposition 2.7. For 3 < k < 6 there are finitely many k-regular graphs 
for which the condition ([23|) is satisfied for an eigenvector of fi\ . 



We next discuss graphs which have an eigenvector satisfying (|23| ) with 
the eigenvalue /i = 2r\ > 2. Recall that by (f25|) r\\k. By counting the edges 
connecting the vertices in two consecutive levels one can show (as for \x = 2) 
that the number of vertices satisfies 

|G| = 2^ r ^n (27) 
Also, since any vertex u\ G G\ has r\ distinct neighbors in Go, 

n > Ji- 
lt is easy to construct examples of regular graphs which have eigenvectors 
with the eigenvalue n > 2 satisfying (p3|); the construction is similar to that 
for u = 2. 

We summarize the previous results: 

Theorem 2.8. Let G be a k-regular graph which has an eigenvector of L 
with an eigenvalue u satisfying (23). Then G is bipartite, n = 21 is an even 
integer dividing 2k, the number of vertices of G is divisible by 2^ k ^ l \ and for 
no > / there exist such graphs with n = 2^ k ^ l ^n,Q vertices. 
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3. Some heuristics for the number of spanning trees 

Here we shall give some heuristic arguments concerning some properties 
that graphs with many spanning trees are expected to have. One such 
property is having large girth; the sequences of cubic graphs described in 
section || seem to support that expectation. We shall give some heuristic 
arguments in its support. 



We start with the following formula from | Mck83 , p. 153] for the number 
r(G) of spanning trees of a fc-regular graph G with n vertices: 

T(G) = , 4(t)(S?+ (^ + _^|Q ( _ iril ,,^ v ', 28) 

where Ik(t) is a function depending on the valency k only (and hence the 
same for all /c-regular graphs), while u r is the number of walks which are 
not totally reducible (i.e. which don't "come" from the infinite tree T^). In 
particular, u r = for r < 7(G). 

The coefficients (-l) r (*) in the sum in ( p8|) are always negative for small 
t, so each nonzero term in the sum decreases the value of the expression in 
brackets in (|2^). Heuristically then, the more u r -s are equal to zero, the 
larger the expression in brackets is in (|28| ) and so the larger is the limit 
which is t(G). Now, if we want to compare t(G) for different graphs in 
those with larger girth 7(G) have more u r -s equal to zero and therefore are 
the natural candidates for large values of t(G), and even more so Moore 
graphs, since they are the unique graphs in having the largest girth. 

Of course, we have only looked at the first nonzero term in the infinite 
series, but we remark that the subsequent u r -s get multiplied by (— l) r (*)fc _r 
which is exponentially decreasing in r. This seems to give at least a partial 
heuristic explanation of the cages appearing in section [|. 

It also seems reasonable to conjecture that graphs with higher connec- 
tivity have more spanning trees. Indeed, one can show that if an edge e is 
a bridge in a one-connected graph G = G\ U {e} U G2 and if G\ (say) is 
two-connected, then one can make an edge switch in G which will increase 
the number of spanning trees (cf. | ]Kel97 ] for the discussion of operations on 



graphs increasing their number of spanning trees). However, it seems to be 
more difficult to prove that for a two-connected graph there exists an edge 
switch which simultaneously increases the number of spanning trees. We 
believe this to be true in essentially all cases, however. 

4. Graphs on a fixed surface 

4.1. The tree number for regular graphs on a surface. In this section 
we shall give an upper bound for the number of spanning trees for 6-regular 
graphs on a fixed surface S which is strictly smaller than the McKay's bound 
(||) for general 6-regular graphs. 

Let G be a /c-regular simple graph on n vertices embedded in an orientable 
surface S g of genus g. Let V, E, F denote the number of vertices, edges and 
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faces of G. By Euler's formula, 

F = E -V + 2-2g > n(k/2 - 1) + 2 - 2g. 

(the inequality becomes an equality if all faces are simply-connected). It fol- 
lows easily that if infinitely many simple /c-regular graphs can be embedded 
in a fixed surface S, then k < 6. We denote by X n ,k,g the set of all simple 
/c-regular graphs embeddable on the orientable surface S g . 

McKay proves in [ Mck83|| that for a random /c-regular graph on n vertices, 



(r(G)) 1/n _> Gk = (k - l^/ikik - 2)) fc / 2 " 1 as n -> oo. We now state the 
main result of this section: 

Theorem 4.1. Let X n f. g be as above. Then 

lim ( sup r(G) 1/n | < a\ < a 6 . (29) 

n ^°° \GeX n , 6 , g J 

To prove the theorem we first establish a proposition estimating the ratio 
of the tree number of G to the tree number of its dual graph G* . We next 
remark that the dual of a 6-regular graph on n vertices is "roughly" a 3- 
regular graph on 2n vertices; an application of McKay's theorem finishes 
the proof. 

Proposition 4.2. Let {Gj} be a sequence of simple graphs on an orientable 
surface S of genus g such that their maximal degrees are bounded. Let t(G) 
and t(G*) be the corresponding numbers of spanning trees, and let nj — > oo 
be the number of vertices of Gj . Then 



-J 29 « ^ « *?■ (30) 



Proof. We first define a mapping from the spanning trees of Gj into (subsets 
of) the spanning trees of G*j as follows: take a spanning tree T of Gj, and 
consider the complement S = Gj — T of T in Gj . It is easy to see ( pgrH ) 
that S is connected, spans G* and that one can get a spanning tree of Gj by 
removing at most 2g edges of 5. We denote by F(T) the set of all spanning 
trees of G* which can be obtained from T by the above procedure. By 
reversing the construction, it is easy to see that the mapping T — > F(T) is 
onto. 

Let |V(Gj)| = rij, and let the maximal degree be bounded above by 
5. Let T be a spanning tree of Gj, and let S = S(T) = Gj — T. Then 
\S\ = \E(Gj)\ - nj + l< nj(5/2 - 1) + 1 and so 



inn <_ fX") « -J 



2;/ 



Adding the inequalities for all T proves one of the bounds in (pOf). The other 
inequality is proved similarly. We remark that the constants in (^) depend 
on 5, k and g only. 

□ 
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Consider now a simple 6-regular graph G (embedded on a surface S g ) 
with n vertices and m = 3n edges, and let G* be the dual graph (of genus 
h < g). G* has m! = 3n edges and at least 2n + 2 — 2h > 2n + 2 — 2g 
vertices; moreover, since G is simple, every face of G has at least three 
sides, so the minimal degree of G* is at least three. It is easy to see that G* 
can be obtained from a (possibly disconnected) graph G\ of maximal degree 
at most 3 with at most 2n + 2 vertices by adding 2h "new" edges between 
the vertices of G%. 



Proof of Theorem |4.l[ . Let G be a graph in X n ^ g . By Proposition 42 



it suffices to estimate r(G*). Consider a spanning tree T of G* containing 
/ of the "new" edges G* — G\, where < I < 2h < 2g. Removing the / 
edges from T produces a spanning forest of G\ with I < 2h components. It 



suffices to establish the estimate (29) for /-component spanning forests of 



G± for each / separately. For I = the estimate follows easily from McKay's 
estimate (||) for cubic graphs on at most 2n + 2 vertices. 

To establish the estimate for I > 1, we first remark that the number of the 
connected components of G\ is at most 2h < 2g; let H\, . . . , H p be denote 
these components, and let Hk have vertices, where q\ + . . . + q p < 2n + 2. 
Each component H = is a graph of maximal degree at most 3. 

Let F be an /-component spanning forest of G\. Let F n Hk = F^ be a 
spanning forest of Hk with Ik components, where l\ + . . . + l p = I. Fj can 
be obtained from a spanning tree of H by removing Ik edges from the tree. 
Each tree has % — 1 edges, so one can get at most spanning forests 

by removing Ik edges from it, so the total number of such forests is at most 



r(H k 



Multiplying the inequalities for all Hk-s we estimate the number of /-component 
spanning trees of G\ by 

ft r(H k ) f[ ( qk ~ ') (31) 

k=l k=l \ Lk ) 

Now, let a sequence Gj G ^6, nj ,g be a sequence of 6-regular graphs on S g , 
and let rij — > oo. Taking the n^-th root in ( [31] ) and letting rij — > oo shows 
that it suffices to prove that 




hmsup t(%) < ^ (32) 



where the subscript j refers to the graph Gj. 

In order to arrive at a contradiction, we shall assume that there exists a 
sequence of graphs Gj on S g such that the limit in ( |32[ ) is strictly greater 
than (T3. Passing to a subsequence if necessary we may assume that the 
number 1 < pj < 2g of components of G±j is actually constant and equal 
to p. We shall order the components Hkj so that qk,j-s are non-decreasing. 
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Passing once again to a subsequence if necessary, we may assume that 
there exists a nonnegative number r < p and a sequence 1 < q\ < qi < 
. . . < q r such that q^j = qu for k < r, and q k j — > oo for k > r. Now 
applying McKay's bound (||) for each k > r (with n = Qk,j) anc ^ estimating 
each of the first r terms in (p2) by an absolute constant leads easily to a 



contradiction, finishing the proof of Theorem 4.1 



□ 

The bound proved in Theorem 4.1 is not optimal. We define the optimal 
bound for 3 < k < 6: 

Jim ( sup (r(G)) 1 /") := a k>g (33) 



n— >oo 



It follows from the proof of Theorem that for g\ < gi 



sup t{G) < n 1 ^ sup t{G) . 

It follows that the limit in ( p3| ) does not depend on g (as long as X n ,k,g 
is not empty for large n, which can happen for small g). Accordingly, we 
define c s (k) to be^| the limit crj. g achieved in ( |33[ ) for some fixed j>0 and 
3 < k < 6. Noga Alon has informed the authors that he can show that 

a s (k) < a k 

(where a k is defined in (0)). 



It follows from the proof of Theorem 4.1 that <r s (6) < a s (3) 2 . One can 
actually show by a similar argument that 

a s (6) = a s (3) 2 

We make a conjecture about the graphs that achieve the limit in ([33]) for 
k = 3,4: 

Conjecture 4.3. The limit a s (3) (respectively, <r s (4)j is achieved by quo- 
tients of the hexagonal lattice ( respectively, the square lattice ) in R 2 . 

It would follow from Conjecture gj (cf. jCDSl , pp. 246-247]) that 



<7 S (4) = e 4G/7r ps 3.21 < 3.375 = ct 4 (34) 

Here G is Catalan's constant. The exponent AG/ir in (j34[) is equal to the 
toplogical entropy 

/ / log(4 — 2cos(27rx) — 2cos(27ry)) dxdy 
Jo Jo 

of the essential spanning forest process on the nearest neighbor graph on Z 2 
( |Solo| ). The integral above is equal to the Mahler measure of the polynomial 
4 — x — 1/x — y — 1/y; it is equal to 

2L'(X,-1) 



3 Here the subscript s stands for surface. 
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where L(x, s) is the Dirichlet L-function associated to the quadratic exten- 
sion Q(V — 3)/Q (cf. |RV| ). It would be interesting to find candidates for 
achieving <r s (5). 

4.2. Jacobians of regular graphs. In this section we want to establish 
an analog for graphs of a result of Buser and Sarnak ( fBSfl ) about the rate 
of growth of the minimal norm of a period matrix of a Riemann surface. 
They have shown that as the genus of a surface increases, the minimal norm 
grows at most logarithmically in the genus. 

Bacher, de la Harpe and Nagnibeda in their work ([BHN, Nag]]) defined 



two lattices associated to any unweighted graph, the lattice of integral flows 
A 1 and the lattice of integral coboundaries N 1 . The quotient (A 1 )*/A 1 is a 
finite abelian group called a Jacobian of a graph (in an analogy to Jacobians 
of Riemann surfaces) . 

We shall now consider fc-regular graphs for some fixed k > 3. Let G 6 
X n the dimension of A}(G) is equal to kn/2 — n + 1. The determinant of 
A X (G) is equal to the tree number t(G), and the (unnormalized) minimal 
norm v(G) of A X (G) is equal to the girth of G. 

We rescale A 1 so that the volume of its fundamental domain is equal to 
1; the (normalized) minimal norm v{G) is given by 

u(G) = 7 (G) r(G) 2/(fen/2 - n+1) . 

From (||) and @ we conclude that for every e > 0, 

6 2/(k/2-l) _ e < r(G)2/(fcn/2 _ n+1) < + £ 

for n large enough. 

We want to study the behavior of v(G) when k is fixed and n — > oo. Using 
the known results about the behavior of the girth, we conclude: 

Proposition 4.4. The minimal norm v{G) of the Jacobian of a k-regular 
graph G € X n ^ satisfies 

v{G) = 0{\og k _ in ) 
This establishes an analogue of the results of Buser and Sarnak. 

5. Some explicit constructions of cubic graphs 



It is known ( pG| , [Kotz| , [Tut]]) that every 3-connected cubic graph may 



be constructed from a tetrahedron by a sequence of edge insertions: given a 
cubic graph on 2n vertices, choose a pair of edges, subdivide them by putting 
a single vertex in the middle, and connect the two subdivision vertices to get 
a cubic graph on 2n + 2 vertices. So, every 3-connected cubic graph G may 
be constructed from the tetrahedron provided we know the "history," i.e. 
the sequence in which insertions were applied to successive pairs of edges. 
Clearly, that "history" is not unique. However, one may hope that if G is 
an extremal graph, among all its histories there exist ones with some nice 
properties. Below we shall give some examples when that is indeed the case. 
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The first example is the following sequence S: starting from a tetrahedron, 
insert edges in such a way that the number of spanning trees is maximized 
after each application. When the number of vertices is small, one may find 
the graphs in S from the table of graphs in [Val]. The 6- vertex graph is 



1^3, the 8- vertex graph is a Mobius ladder (an 8-cycle with the pairs of 
''opposite" vertices connected, cf. [ Big93| , §3e]); it has 392 trees. The 10- 



vertex graph is the Petersen graph P (2000 trees), and the 12- vertex graph 
is obtained from P by applying A to a pair of edges at distance 2 from each 
other (this is graph # 84 in JVa]||; it has 9800 spanning trees). The last 



graph also appears in [Hoi]. From the table of graphs in | |Val| one can check 
that the graphs in the sequence have the most trees among the cubic graphs 
with the same number of vertices. 

The rest of the graphs in the sequence S are shown in the Appendix; they 
were obtained by a computer program. It is remarkable that the sequence 
contains all (3,g)-cages for 3 < g < 8. Namely, the 14-vertex graph in S is 
the Heawood graph, the 24- vertex one is the McGee graph, and the 30- vertex 



one is the Tutte-Coxeter graph (cf. [ Won ] or | HS , Ch. 6]). Unfortunately. 



the graphs in S with 32 or more vertices seem to be less interesting. The first 
few graphs in the sequence appear in |AS[| (but without the reference to the 
tree number); however, the inclusion of the 12- vertex graph in the sequence 
and the graphs with 16 or more vertices appear to be new. The authors hope 
that similar constructions may lead to discovery of new extremal graphs .[] 

Edge insertion has its disadvantages: for example, certain cages of even 
girth have the property that every two edges are contained in a common 
shortest cycle ( |KT| ); inserting an edge will then certainly decrease the girth. 
So, one might try adding several new vertices and edges simultaneously. 
For cubic graphs, one can try subdividing several edges of an "old" graph 
(with one, two or more points on an edge) to get an even number of "new" 
vertices of valence 2, then choosing a perfect matching for the new vertices 
and connecting pairs of these vertices according to the perfect matching to 
get a new cubic graph. 

Several cubic cages can be obtained by taking uniform subdivisions of 
smaller graphs and inserting edges 

• The Petersen graph can be constructed from a tetrahedron by subdi- 
viding its 6 edges; 

• The Heawood graph can be obtained from the unique loopless two- 
vertex cubic graph (G2) by putting 4 new vertices on each of its 3 
edges; 



6 



Other sequences of cub ic grap hs might be constructed using other graph invariants 



(e.g. mean distance - cf. [CCMS| - instead of the tree number) in the algorithm; one 
can also start from a "large" cubic graph and delete an edge together with its endpoints 
so that a given invariant is maximized. For regular graphs of degree k > 4 one can use 
operations similar to the edge insertion to increase the number of vertices and construct 
sequences of fc-regular graphs b y ma ximizing a certain invariant. 



7 Similar examples appear in [AS] 
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• The Tutte graph can be obtained from 2K3 3 by subdividing each of 
its 18 edges. 

Appendix: a sequence of cubic graphs; constructing cubic 

CAGES. 

This section describes a sequence of cubic graphs generated (by a com- 
puter) from the tetrahedron by edge insertion where at each stage the in- 
sertion was chosen so that the number of spanning trees of he new graph 
was the maximal possible. The first six graphs are shown in Figure 1; the 
endpoints of the new edge inserted into a cubic graph on 2n vertices were 
numbered 2n + 1 and 2n + 2. 
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Figure 2. The 18-vertex graph and the McGee graph. 
The Tutte-Coxeter graph (as it appears in S) is shown on Figure 3. 
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Figure 3. The (3, 8)-cage. 

In Figure 4 we show how to construct Petersen and Heawood graphs by 
edge insertions into the tetrahedron and the graph G2 respectively. The 
vertices of K4 and G2 are shown in black. 





Figure 4. Petersen and Heawood graphs. 

In Figure 5 we show how to construct the Tutte graph by edge insertions 
into 2X3,3. The vertices and the edges of 2^3,3 are shown in black; the 
vertices of the two copies of ^3,3 are marked A and B respectively. 
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Figure 5. The Tutte graph obtained from 2K^^. 
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